The Fermi edge singularity in the SU(N) Wolff model 



Balazs DoreQ 

Department of Physics, Budapest University of Technology and Economics, H-1521 Budapest, Hungary 

(Dated: February 2, 2008) 

The low temperature properties of the SU(N) Wolff impurity model are studied via Abelian 
bosonization. The path integral treatment of the problem allows for an exact evaluation of low 
temperature properties of the model. The single particle Green's function enhances due to the 
presence of local correlation. The basic correlation function such as the charge or spin correlator 
are also influenced by the presence of impurity, and show local Fermi liquid behaviour. The X-ray 
absorption is affected by the presence of local Hubbard interaction. The exponent is decreased 
(increased) for repulsive (attractive) interactions. 
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I. INTRODUCTION 



Local interactions have played an important role in explaining the behaviour of localized impurities, and in under- 
standing the effect of electron correlation in lattice theories^. The most studied magnetic impurity models are the 
Anderson and Kondo model, including their various generalizations 2 . Their local Fermi or non-Fermi liquid behaviour 
in the low temperature regime is inferred from Bethe ansatz method, but dynamical properties are still very difficult 
to be calculated. Beyond the reliable renormalization group calculations^, one is left with bosonization to gain insight 
into the low energy properties of such models. 

One of the simplest description of a single impurity in a non-magnetic metal is based on the Wolff model, where 
electron correlation is still presents^^. Also it is the basis of studying the effect of Coulomb interaction on resonant 
tunneling through a single quantum levelZ^. The model consists of conduction electrons interacting with each other 
solely at the origin through the Hubbard interaction. It is related to the Hubbard model in the same way as the single 
impurity Anderson model is related to its periodic version. 

Several extensions of this model are possible including extended interactions or by changing the properties of 
conduction electrons. In the present work we study a most natural generalization of the usual Wolff model: we 
consider fermions with SU(N) spin index or alternatively we take the ensemble of spin and orbital degrees of freedom 
into account by the N index instead of the usual two-valued spin index. The additional degrees of freedom can be 
realized through orbital degeneracy, for example, as in Mn oxides^. As a result, the additional degrees of freedom 
can be called flavour or color index. In similar models, the inclusion of orbital or flavour index led to non-Fermi 
liquid behaviour, such as the two channel Kondo modelifi*ii, or enriched our picture concerning the Mott-Hubbard 
transition as in the SU(N) Hubbard modeli 2 *^. 

The purpose of the present investigation is to see, how the interaction between N different species of electrons at 
a single site influence the low energy properties of the system. In the SU(N) Hubbard model with one particle per 
site, the Mott phase is reached at a finite value of the Hubbard interaction^ for N>2, as opposed to the zero critical 
interaction for the SU(2) casei. In the SU(2) Wolff model, the crossover to strong coupling region is predicted to be 
symmetric with respect to the sign of U£*&. Here we shall see how this condition is modified for SU(N) symmetry. 
The use of Abelian bosonization allows for an exact evaluation of the basic correlation functions, which show the 
instability toward local magnetic moment formation (i.e. Kondo regime), as the interaction parameter is increased. 
The critical interaction for repulsive U is independent of the value of N, but for attractive interaction its value is 
suppressed by 1/(N-1). The low energy behaviour of the system is characterized by a local Fermi liquid, obeying the 
well-known Fermi liquid relations (such as the Korringa relations) of the Anderson impurity model. 

The X-ray response of the SU(N) Wolff model is evaluated with the use of boundary condition changing operators^ 
and is qualitatively similar to a system with correlated electrons over the whole lattice (namely the Hubbard model)i^. 
The exponent characterizing the scaling dimension of the deep level electrons decreases with increasing Hubbard 
interaction. 
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II. FORMALISM, MODEL HAMILTONIAN 



The Hamiltonian describing N different species of electrons interacting only at the origin is given by: 



N 



N 

J dx*+(x)d x * m (x) + E : *+ (0)* m (0) : +^ £ : ¥+(0)* m (0) :: *+(0)*„(0) 



n— l,n^m 



, (i) 



and only the radial motion of the particles is accounted for by chiral (right moving) fermion fields^. The model can 
be bosonized viai*i& 



and the bosonized Hamiltonian takes the form as 



H = v f d I ^(9 I $ m ( I )) 2 + ^4$„ i (0) + ^ d x $ m (0)d x $ n (0). 

_ OQ m V m n,m=l,n=£m 



(2) 



(3) 



By introducing charge and spin fields as 



v m— 1 

*n,.(aO = . = [ E - n$„ +1 (.T) ] , 



(4) 



(5) 



ri = 1. . . N — 1, the Hamiltonian separates into different sectors: the spin sector is described by N — 1 identical 
decoupled bosonic modes, and the charge sector transforms into a similar massless bosonic mode as 



N-l 



dx(d x $ m , s (x)) 2 - ^-(d x <P m , s (0)) 2 

Z7T 



//,. = <• / dx(d x <P c (x)) 2 + E ] /^d x <s> c (o)+ {N 2 ^ )u (d x <s> c (o)) 2 - 



(6) 



(7) 



Since the Hamiltonians describe scattering of massless bosons on a point-like "impurity", the bosonic correlation 
functions can readily be calculated. Hence in the followings it suffices to study the single Hamiltonian 



H = v dx{d x $>{x)f + Ed x $(0) + U(d x <P(0))< 



(8) 



describing both sectors with the appropriate parameters. When evaluating the generating functional in terms of path 
integral, one can benefit from the fact that the Hamiltonian depends only on d x <5>(x), and not separately on the dual 
fields 4>{x) and 9(x). As a result, the path integral can be formulated after identifying the "space coordinate" with 
$(x) and the "momentum" conjugate to ifcll with —d x $>(x) from the commutation relation of the $ field with itself 
([$(&), = isgn(x — y)/4). The obtained action is for chiral bosons, similar to the one investigated during the 

study of edge statesA. Then the generating functional is obtained as 



with 



S = dr' 



-i I d T ,<Z>(x,T')d x <Z>(x,T')-v I dx(d x $(x,T')) 2 - Ed x $(0,T') - z(t') 2 -2z(T , )V z Ud x <$>(Q 1 T') 

+tp($(0,r)-*(0,0)), 



(9) 
+ 

(10) 
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where the interaction term was decoupled by a time dependent Hubbard Stratonovich field z(t). For /i = 0, it gives 
the partition function, for /x = y/An/N it equals to the local Green's function at the impurity site (G(/x) = Z(p)/Z(0)) 
coming from either the charge or one of the spin sectors. After Fourier transforming the fields as 

$(x, t) = ^ exp(-iw m r - iqx)$(m, q), (11) 

m , q 

z(t) = ^ exp(-iw m T)z(m), (12) 
m 

Lu m = 2m7rT is the bosonic Matsubara frequency, the action is transformed to 

S = ^2(-iqcj m - vq 2 )$(-p)$(p) + iEy^g$(0,g) - ^ z(-m)z(m) + 2fy/ z J7^ z(-m)q$(p)+ 

p q m p 

+ i/x(exp(-iw m r) - l)$(p), (13) 
p 

where p = {m, g}. The $ field can be integrated out, and the effective action reads as 

1 

2 



= - J2 <-m)z{m) + - ^ vi-p)Gv{ P ), (14) 



m p 

where 

G _1 = vq 2 + iqu m , (15) 
r](p) = iEqSmfi + 2iy/—Uz(—n)q + i/i(exp(— iu> m r) — 1). (16) 



Integrating over the z field gives for the free local Green's function (where one has to substitute /i = y 4ir /AT) in the 
large r limit 

2 T 1 _ cos(w m T) \ / T y 2/A * 



G(r,0)~exp — > — " ~ rj^rr • (17) 

y 2 p V Q +i<l UJ m J \sm{7rTT) J 



As to the f7 dependent part, it is obtained as 



in the large r limit with 

ioj + vq 



7W = E5mI^. ( 19 ) 
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AK0 = ^-^-, (20) 



Uq = irv/2no, uq is the average density per spin channel in the homogeneous case. For \x = 0, the change of the free 
energy of Eq. © in the presence of U is evaluated as 



F = T EMi + ^ m ))-i I f^ yi (21) 



which gives for the total specific heat 



C(T) = — ^ . (22) 



The presence of local Hubbard U enhances the Sommerfeld coefficient for attractive interaction, and decreases it for 
repulsive one. 

Another way to derive the same results would have been to follow the diagrammatic approach. Due to the presence 
of a single scattering center, the determinations and solution of the corresponding Dyson equation would have led to 
the same exact expressions for the Green's function and free energy. 
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III. CORRELATION FUNCTIONS OF THE SU(N) WOLFF MODEL 



By inserting the appropriate parameters into these results, we are able to determine the low temperature-low energy 
properties of Eq. Q • First of all, the specific heat including the impurity contribution reads as 

*T ( U (N-1)U \ 

as opposed to its normal state value C(T) = NttT/Qv. In Fig. we show the Sommerfeld coefficient, Cw(T)/C(T) 
as a function of U. The main quantity of interest is the local electron propagator. Putting the above results together, 
one finds for the single particle electron Green's function at the impurity site in the 0<t« 1/T limit 

G m[ r) = -<* m M)*+(0,0)> = - ^^)W^W^SW^W^ ' (24) 

where U = irv/2n is the critical Hubbard interaction in the SU(2) Wolff impurity model, n is the average density 
in the homogeneous system. For general N>2, however, the critical value of U is suppressed by 1/N as seen from the 
denominator of Eq. 1241 When U exceeds this value, we expect local magnetic moment formation, and bosonization 
breaks down^&. The one-body potential, E does not appear in Eq. Q24fl to leading order in r, only renormalizes 
the t -2 correction. The r — > oo limit determines the low energy behaviour of the local density of states. It remains 
constant, i.e. Fermi liquid like, but enhances due to the presence of local Hubbard U, as can be seen in Fig. ^ 
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FIG. 1: The Sommerfeld coefficient (left panel) and the density of states at the Fermi energy (right panel) are plotted as a 
function of local Hubbard U for N = 2, 3, 4 and 8 from edge to center. As is seen, the critical value of the interaction depends 
on its sign, namely U c = Uo for repulsion and — U c — Uo/(N — 1) for attraction. 



In the SU(2) case, the Sommerfeld coefficient and the local residual density of states have the same form, because 
the Hamiltonian splits into two sectors somehow corresponding to the original two species of electrons. For N>2, 
however, we have N different species, but only have two distinct sectors with different weights: N-l for the spin sector 
as opposed to the single charge sector, resulting in different behaviour in the specific heat and Green's function. 

The other, r — > + limit determines the particle density as 

\/NE 

The computation of the local density-density correlation function leads to 

1 + (N - l)(7xoM 

where Xo( w ) is the local density correlation function in the homogeneous case, Xo(0) = Due to Ward identity^, 

the exact susceptibility is identical to the RPA result in the SU(N) case as well&. The above result can further 
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be verified by the fluctuation-dissipation theorem. Namely, the charge response function is related to the charge 
correlator by (n 2 ) — (n) 2 = TXcharge- Here n is the charge density, and its momenta can be obtained by deriving 
the partition function with respect to E. At low frequency, where bosonization is expected to be valid 1 , it is well 
approximated by 

N . u NU 2 

Xcharge^) « Uq + ( N _ 1)f/ + ' 2 7TV 2 (U + (N - 1)U) 2 ' ( ' 

Its frequency dependence is shown in Fig. [2 for N=4 and for various values of U. As is seen, the low frequency part 
is influenced more significantly for attractive interactions. As U decreases below zero, its imaginary part becomes a 
Dirac delta function at the origin, signaling the breakdown of bosonization. Similarly, the susceptibility in the spin 
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FIG. 2: The imaginary part of the density-density response function (left panel) and the imaginary part of the spin susceptibility 
are shown as a function of frequency for N=4, U = 0.9Uo, 0.5(7o, 0.25(7o, (dashed line), — 0.2£/o and — 0.33£/o with increasing 
(decreasing) initial slope at uj = for the charge (spin) response. W is the bandwidth assuming a sharp momentum space 
cut-off. 



sector can be evaluated as 

(jv-i)xoH 

Xspin{U) = j-j — f-^-, 28) 

1 - Uxo(u) 

which again reduces to the RPA results. This is N-l times the spin susceptibility of the SU(2) Wolff model, stemming 
from the N-l massless bosonic modes in the spin sector. The low energy spin response is evaluated as 

N-l uj (N-1)U 2 
Xspin M « ^— jj + ^ {Uq _ u? ■ (29) 

Compared to Eq. (I27II . for repulsive interaction the spin response increases faster in energy than the charge one due 
to the denominator Eq. H29|) . as can be checked directly in Fig. [21 Similar results were found in the SU(2) case&. 
Exactly at the critical U deduced from the Green's function, the charge or spin susceptibility diverges, signaling a 
possible transition or crossover to strong coupling case. 

IV. X-RAY ABSORPTION 

The basic Hamiltonian describing our previously studied system and deep level electrons (d) is given by 

H deep = H + E Q J2d+d n + V *m(0)*m(0)#<k, (30) 
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and the coupling of the metal to the X-ray field, describing the transfer of a deep level electron to the conduction 
band is given by 

E x = WJ2 *+(0)d m e- iart + h.c. (31) 

m 

The quantity of interest is the X-ray scattering rate, which is written in the form of the Fermi golden rule as 

I{lu) = 2\W\ 2 lmS{uj), (32) 

S(t) = (d+(t)* m (t, 0)^(0,0)^(0)). (33) 

As is well known, the number of d electron is conserved ([d + d, Hdeep]), it is not a truly dynamic quantity, so one can 
reformulate the problem in terms of time dependent core-hole potential. After bosonization, the core level electron 
couples solely to the charge sector as 



H deep = H + E a d+d m + vJ^^2d+d m d x <f> ( 



(0). (34) 



From now on we focus on the excitation of a single d m electron, and forget about the presence of the other deep 
level particles. They can only renormalize the threshold frequency but not the exponent what we are interested 
Following the strategy of Schotte and Schottei^, we introduce two Hamiltonians, one when there is no hole 
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(Hi = H) and another, when conduction electrons feel a scattering potential (Hp = H + Eo + V^2^ l=1 ^+ (0)^ m (0)). 
Similarly to the case of a normal Fermi liquid, there exist a unitary transformation (the boundary condition changing 
operator—) relating the two Hamiltonians as 

H F = U + HiU, (35) 
[7 = cxp ( z^l$ c (0) ) , (36) 



where 



VT ^ (37) 



2v U + (N-1)U 
With the use of these, the d electron Green's function can be calculated from 

D(t)=-{d m (t)d+) = (U+(t)U). (38) 

Its time dependence can easily be determined from the asymptotics of the correlator of the charge sector using Eqs. 
DIlandEaas 

D(t) ~ i-W*?. (39) 

As is seen, the time decay implies non- Fermi liquid behaviour, characteristic to the X-ray problemM. The Fourier 
transform of the deep level electron can be obtained, from which the d-electron density of states follows as 

, . sin(7r(l - (5/n) 2 ) , , 

The original Dirac delta density of states situated at Eq transforms into a power-law divergence with a finite tail. 
Compared to that of a normal metal without Wolff impurity, the exponent 5 decreases for large number of channels, 
while in a normal metal the exponent grows with y/~N. Also at fixed TV, repulsive (U > 0) interaction suppresses 
5, while in the attractive case the exponent increases rapidly. It is worth mentioning, that since only the charge 
channel is involved in the X-ray Hamiltonian, the singularity present in the Green's function (Eq. 1241) is absent here 
for repulsive interactions. The two-particle Green's function, S(t) can be evaluated similarly to D(t), namely 

S(t) = (U(t)9 m {t,0)9+(0,0)U+). (41) 

Using the asymptotic formula for the correlator of the $ field, it reads as 

S(t) ~ t- aiSM'+N-l)/N_ (42) 

Since only the charge degree of freedom is involved in the conduction-deep level electron interaction, the scattering 
phase shift comes with an additional factor of 1/N, overwhelmed by the spin degrees of freedom with zero phase shift. 

Similar phenomenon occurs in the X-ray edge absorption, when we extend the local correlation of our model to 
the whole lattice, namely in the SU(2) Hubbard model. The X-ray exponent was found to decrease for repulsive 
interaction^, and increase for attractive one. We believe, that the qualitative dependence of local quantities on 
electron correlation is well described by the single impurity Wolff model. 
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V. CONCLUSIONS 



We have studied the low energy properties of the SU(N) Wolff impurity model via Abclian bosonization and the 
path integral treatment of the generating functional. The system decouples into N-f identical spin sectors and one 
single charge sector, involving only chiral bosonic fields. The single particle Green's function possesses Fermi liquid 
type asymptotics, but its weight can strongly be modified by the local Hubbard interaction. The spin and charge 
sectors become inequivalent, signaling a transition to the strong coupling phase at —Uq/(N— 1) for attraction and at 
Uq for repulsion. For small interaction, the system exhibits local Fermi liquid behavioui*2> in all evaluated quantities. 

When interaction with deep level electrons is considered, the X-ray response of the model can be explored. It can 
be formulated similarly to the original work of Schotte and Schott ob 19 , exploiting the fact, that the number of deep 
level electrons is a conserved quantity. The deep level electron Green's function can be obtained using the boundary 
condition changing operators^, reflecting the influence of both the local Hubbard interaction and the presence of N 
species of electrons. The X-ray exponent is decreased (increased) for repulsive (attractive) Hubbard interaction. 
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